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Abstract 

We apply the AdS/CFT holography to the simplest possible eternal time machine solution in 
AdS3 based on two conical defects moving around their center of mass along a circular orbit. Closed 
timelike curves in this space-time extend all the way to the boundary of AdS3, violating causality 
of the boundary field theory. By use of the geodesic approximation we address the “grandfather 
paradox” in the dual 1 + 1 dimensional field theory and calculate the two-point retarded Green 
function. It has a non-trivial analytical structure both at negative and positive times, providing 
us with an intuition on how an interacting quantum field could behave once causality is broken. 
In contrast with the previous considerations our calculations reveal the possibility of a consistent 
and controllable evolution of a quantum system without any need to impose additional consistency 


constraints. 
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I. INTRODUCTION 


Solutions to the equations of General Relativity that describe space-times containing 
closed timelike curves (CTC) have attracted significant interest as they revealed at least 
hypothetical theoretical possibility of travelling in time. Since the renowned publication by 
Kurt Godel [1] a number of causality violating solutions in GR as well as in modified theories 
of gravity have been constructed, among which we can name the Tipler-Van Stockum time 
machine generated by axially rotating distribution of particles [2], [8], the Morris-Thorne- 
Yurtsever transversable wormhole [4] 5], the Gott time machine based on moving conical 
defects |6], the Ori dust solution [7], and the solutions in f(R) theories of gravity [8] and 
theories with non-minimal matter-curvature coupling [9]. 

All questions about physics of time machines that could be posed in principle fall into 


three general categories: 
e Is there a physical way to create a time machine? 
e Is there any time machine solution that can be stable? 


e What dynamical behaviour would a physical system experience evolving in a time 


machine background? 


None of the questions have yet received a definite answer. 

The answer to the first question is believed to be negative. Extensive analysis of particular 
time machine solutions has demonstrated that in order to create a space-time with CTC 
one needs matter that violates strong, weak or null energy conditions of General Relativity 
(different solutions require violation of different energy conditions), and only eternal time 
machines can exist [10|{I1]. However we can not be sure that all matter in the Universe obeys 
these conditions. For instance, there are a number of models of the dark energy violating 
the null energy condition 15], and this provides a way to by-pass the no-go statement. 

The second question was raised by Hawking in [16], where he conjectured that a space- 
time with CTC can be stable only on classical level, but will be unavoidably destroyed by 
quantum fluctuations of the metric. The real universal proof or refutal of the conjecture 
can be obtained only within a framework of a complete theory of quantum gravity. String 
theory opened a possibility to check the chronology protection condition in specific cases. 


In authors have shown that appearance of closed timelike curves in a certain (O-plane) 
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orbifold background would cause a Hagedorn transition that restructures the space-time 
transforming it into a chronologically safe configuration. So this result can be considered as 
a very accurate and nice supporting evidence in favour of the Hawking conjecture. On the 
other hand in [18H20], it was demonstrated that the Gédel type solutions can be smoothly 
embedded in the context of string theory. Closed timelike curves in that case are hidden 
behind the so called holographic screens and do not violate causality in the rest of the 
space-time. Thus the chronology is protected, but structure of the CTC remains unbroken 
by quantum effects. An intriguing observation has been made by authors of [21] and [22], 
that from the point of view of the AdS/CFT correspondence, the existence of CTC in the 
bulk can be related to negative or exceeding one fermionic probability in the boundary field 
theory. 

The third category includes various types of the “grandfather paradox”. For a clas- 
sical wave equation on a non-globally hyperbolic space with CTC the possibility of self- 
consistent dynamics was demonstrated in [23] [24]. Classical mechanical billiards and their 
self-consistency conditions have been studied in [25]. A basis of states of a free quantum 
field theory in the Gott time machine has been constructed in , where it was shown that 
the causality violation leads to an emergence of an effective non-unitary interaction in the 
theory. Non-unitarity of interacting field theories in time machines was analyzed in [27]. 
Some authors even argued that evolution of a physical system along closed timelike curves 
can be studied experimentally by mean of simulation of emergent gravity in metamaterials 
28], or a qubit interacting with an older version of itself [29]. 

However the question about properties of an interacting quantum field theory in a time 
machine background remains open, though the real “grandfather paradox” can take place 
only in a self-interacting system. When the notions of time ordering and unitarity are 
absent from the very beginning, it is unclear how to formulate an interacting field theory. In 
this paper we address this problem and by use of the AdS/CFT correspondence provide a 
constructive solution to it. Here we stand on the position that even if the presence of CTC 
causes breaking of unitarity in the boundary field theory [22], it should not be regarded as a 
big problem as long as we can formulate a prescription how to solve the theory. When one 
is trying to get an insight into physics of paradoxical systems, it is not very useful to rely 
on the “common sense” intuition and corresponding fundamental principles. 


The AdS/CFT correspondence provides an elegant way to address the paradox. In the 
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large N limit it relates quantum field theories to a classical gravity, and thus we can study 
properties of a quantum theory in the CTC background just by careful analysis of the dual 
Riemannian geometry, without any need to formulate special quantization rules that would 
be valid in the case of broken causality. 

The paper is organized as follows. In the next section we introduce a simple set up 
for the time machine in AdS3. In Sec. [III] we discuss geodesic structure of the spacetime, 
and suggest that it could lead to non-trivial effects in the boundary field theory. In Sec. 
we introduce a notion of timelike quasigeodesics that will be then used for connecting 
timelike separated boundary points. Finally, in Sec. [V] we formulate a precise algorithm for 
the Green function evaluation, provide the results of numerical simulations, and discuss the 


related phenomenology. 


Il. TIME MACHINE IN AdS3 


The eternal time machine solution in AdS has been suggested by Gott and DeDeo in 
(for similar solutions containing CTC but collapsing into a BTZ black hole see ). Here 
we briefly recall its structure closely following the original text. 


The three dimensional global anti-de Sitter space-time can be thought of as a hypersurface 


—X¢-—XZ4+ X74 XZ =-1, (1) 











embedded in a four-dimensional flat R?? space-time with a metric: 





ds* = —dX@ — dX3 + dX? + dX3. (2) 
In the Schwarzschild coordinates the embedding formulas are 


Xo = V1+ R?2 cost, (3) 
X3 =V1+4 R?sint, 
X,= Reos¢, 
X2=Rsind, 
where R € (0,00), € (—00, 00), ¢ € [0, 27). 


The induced metric is then 


dR? 


ist Rd¢”. (4) 


ds* = —(1+ R’)dt? + 




























































































FIG. 1. Two equivalent unfoldings of the AdS3 spacetime with a conical defect a = 3m. The 
larger part of the spacetime outside of the faces is to be cut out, and the faces are identified. 
The only physical space is the narrow region between the faces. To construct a time machine the 


twisted unfolding is more convenient to use. 


A massive particle put into a three dimensional space-time removes a wedge with an angle 
deficit proportional to the mass of the particle, and edges (faces) emerging from this point- 
like particle. Points on the opposite edges of the wedge are identified, and the resulting 
space-time contains a conical defect, Figfi{a). When we are looking at the unfolding of 
the conical defect, coordinate locations of the edges do not have an independent physical 
meaning, and we are free to rotate them preserving the angular deficit. For our purposes it 
will be convenient to make the cut out “pizza slice” twist in time with a constant angular 


velocity in the reference frame of the massive particle, making a full rotation in a period 27, 
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Fig/i{b). Then for the trailing and leading faces of the wedge in the embedding coordinates 
we get: 

Xi = V14+F2cost X}=VJ1+ BR cost, 

Xi= V1+Rsint X$=/1+4+ Rsint, 

Xt = Reos(t —a/2) X! = Reos(t+a/2), 

Xf = Rsin(t—a/2) X$ = Rsin(t+a/2). 


Here a is the angular deficit of the conical spacetime. Integrating the spacetime stress-energy 
tensor over the angle, we can deduce that the effective mass concentrated in the interior of 


the bulk is 


1 a 
oe "3G ' Tere’ (5) 





Here the second term is the mass of the point-like source, and the first one is the contribution 
from the negative AdS curvature. 

For a single static conical defect we can not make its angular deficit a larger than 27. If 
the mass of the point-like particle exceeds the limit a = 27 (M > 0), the resulting space-time 
will be rather a BTZ black hole instead of a naked conical singularity [82]. 

However, if we boost the massive source, the removed wedge is effectively getting 
“squeezed” from the point of view of an external observer at rest (see Fig. [2{a)). This 
provides a room for a second conical defect with a deficit angle 6 such that a+ 6 > 2m. 
In other words, relative motion can support the system of conical defects, preventing them 
from collapsing into a black hole. This will be the essence of the DeDeo-Gott construction. 

Consider a system of two identical conical defects undergone two opposite Lorentz trans- 


formations, in the (Xo, X1) and (X3, X2) planes of the embedding space each: 


coshw 0 sinhy 0 
0 cosh O sinh 

Ay = Aj} = . : (6) 
sinhw 0 coshy 0 


0 sinhw O - coshw 


In the three-dimensional coordinates of the AdS spacetime these Lorentzian transformations 
correspond to SO(2,2) isometry transformations. 
It can be shown that in the coordinates of global AdS3 these defects move along the same 


circular orbit R = const with a constant velocity, always being at the opposite points of the 
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FIG. 2. (a) A boosted conical defect in AdS3. Faces of the wedge are deformed, and the identifi- 
cation occurs between points with different time coordinates (in the centre of mass frame). Here 


a = V3, i = 1. (b) The DeDeo-Gott time machine. 


orbit, Fig{2{b). A conical defect sits at R = 0 in its rest frame, ie. 
Xo = cost, X; =0, X» =0, X3 =sint. (7) 
In the boosted frame 
Xo = coshy cost, X; =sinhy cost, X2 = sinhwsint, X3 = coshwsint, (8) 


hence 


R= 4/X?4 X32 = |sinhyI, (9) 


Xo ~~ sinhwcost 
= = sd o 1 
PR jsinh | sgnw-cost (10) 





cos @ = 
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In the rest frame of a wedge, the points on its edges are identified at equal coordinate 
times. However if we boost it, from the point of view of an external observer this identi- 
fication would occur at different times leading to time jumps for a particle moving around 
the conical singularity. In the case of a single conical defect its boost can be regarded as a 
global coordinate transformation of the space-time, which obviously can not cause any new 
physical effects. However, with two defects moving relative to each other, the relative time 
jumps become a physical effect that can not be eliminated by a (proper) choice of coordinate 
system. These time jumps allow for the existence of CTC. 

Existence of closed timelike curves in this space-time can be demonstrated by looking at 
the identification of the edges near the boundary of AdS3 (at R — oo). We refer the reader 
to for a detailed discussion, here we just quote the result. Speaking in terms of the 
unfolding of the two-conical space-time, when a timelike particle living on the boundary of 


the AdS cylinder hits an edge of one of the two wedges, it undergoes a time and an angle 








jump: 
_ sin(a/2) tanh wy 
At = 2arctan (, poco 7 : (11) 
7 sin(a/2) 
Ad = 2arctan (5 rs 2) ’ (12) 


It can be shown that At+ Ad =a. If a > 7, the world line of the particle becomes a closed 
timelike curve, and thus the space-time is a time-machine, see Fig [3] 

In the regime when this limit is not exceeded, and CTC are not present, the space-time 
has been studied in detail from holographic point of view in [33], but the case of broken 
causality has not been addressed. 

In the next sections we will study geodesic structure of this time machine and explicitly 
show that in presence of the two orbiting conical defects we deal with a highly-nontrivial 
lensing of geodesics, and this reflects on the structure of two point Green’s functions of the 


dual boundary theory. 


Ill. ENTWINEMENT OF GEODESICS AND CAUSALITY VIOLATIONS 


When conformal dimension A of a boundary operator in AdS/CFT is very high, the 


corresponding two-point Green function can be derived in the geodesic approximation [32]: 
G(A, B) = e448 , (13) 
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FIG. 3. A timelike particle moving along the boundary experiences a jump in time and angle when 
it hits a face of one of the wedges. If At + Ad > 7a, closed timelike curves become possible. The 
gray strips are the cut out part of the boundary, and the white strips are the physical part of the 


boundary. 


where Lap is the length of a geodesic connecting boundary points A and B. If there are 
more than one geodesic between A and B, they can give additional contributions to the 
propagator. This is the case for the DeDeo-Gott time machine geometry, and here we 
address possible outcome of this in details. 

Consider two arbitrary points A (t,, @;) and B (t2, ¢2) located in the physical (unremoved) 
part of the AdS3 boundary. Having two rotating conical defects in the bulk makes the 
structure of possible geodesics connecting!] A and B very nontrivial, so we should find a way 
to calculate their contributions to the two-pont Green’s function G(A, B). Let us shoot a 
geodesic from the boundary point A to the point B. Before it hits the point B it can undergo 
a number of “refractions” on the faces of wedges, winding around either of two conical defects 
clockwise (if it hits the leading face of the wedge head on) or counterclockwise (if it overtakes 


the trailing face of the wedge from behind). For example, schematically a typical geodesic 





' When A and B are timelike separated we encounter some subtleties caused by the fact that in the AdS 
space-time a timelike geodesic can not reach the boundary. These issues will be commented further on, 


but the general point of view described in this section remains unchanged. 
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may have a structure (see also Fig. 
A>+W, ~Wi > WwW; > Wi, > B, (14) 


where Writ stands for the act of clockwise/counterclockwise winding around the 1-st or the 
2-nd wedge respectively. 


So, formally the Green’s function in the geodesic approximation is given by 


G(A, B) = Fy Meets Aa (15) 
n=0 {W1---Wn} 


where the second sum is taken over all different entwinement structures corresponding to 
the same number of windings, and the first sum is taken over all winding numberg?} 
It is easy to see that for a given number of windings N the maximal possible number of 


topologically different geodesics is 
a oe (16) 


The first winding act can be of four different types. But for each of the next steps, if a 
geodesic wrapped around a conical defect, for example, clockwise, then on the next step it can 
not go in the opposite direction and wrap around the same conical defect counterclockwise. 
It means that in the sequence of windings the winding act W; can be followed (at least 
hypothetically) by W;, W7;, W7;, but not by W,. 

As we will see further, for a given pair of boundary points (A,B) not all sequences of 
entwinements are physically realistic and can contribute to the sum (15). 

To find a proper prescription for the lengths of the non-trivial winding geodesics let us 
discuss in detail a particular example. 

Take a look at Fig [4 The length of the presented geodesic is a sum of lengths of its’ 


composing arcs: 
L(A + W, WW; Wi, > B) = Lao, + £Lo.03 + £o.0; + Logo, + Logs. (17) 


We will refer to points C; as complementary points, and O; as refraction points. Coordinates 
of the refraction points O; are to be found from coordinates of A and B, and that can be 


easily done step by step. 





? This idea of entwinements in holography has been introduced in [34], but there it was related to a concept 


of entanglement entropy “shadows” rather than to subleading contributions to the propagator 
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FIG. 4. Schematic constant time projection of a typical geodesic connecting points A and B in 
the time machine. Red curves are for leading faces of the rotating wedges, and blue curves - for 
trailing faces. The entwinement configuration in this particular case is W, WiWiWi according 
to the notations introduced in the main text. All shown points in principle can have different 
time coordinates (here we schematically project them down to a single time section, so the curves 
the geodesic is made from should be taken only as an approximate artistic representation). On 
picture (a) the acts of entwinement and identifications are shown explicitly. A and B belong to 
the physical unremoved part of the spacetime, and O; are the points where the geodesic undergoes 
“refraction” on the wedges. Picture (b) demonstrates the idea of complementary points C; located 


in the removed part of the spacetime. 


Focus on the first refraction on the wedge, W; . The points of interest are {A, O1, C1, C2, Oo, C3}. 
Here we neglect for a while the second wedge, so we do not consider the point O3 at all, 
and we treat C3 as a physical point (not just as a point in the complementary “removed” 
space), see Fig {a). These six points can be regarded as a result of boost transformation 


A; applied to the wedge. We can “unboost” the wedge and find static pre-images of these 
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FIG. 5. The left picture demonstrates W, refraction of a geodesic on the first conical defect. The 
right picture demonstrates how it looks like in a reference frame of the conical defect. Again, points 


generically belong to different time slices, and the representation is purely schematic. 


points (see Fig 5{b)). For the “unboosted” points the following relations trivially hold 


Cy = Aj1C, = Aj1A — (0,a) = A-(0,a), (18) 
C3 = Ay'C3 = Aj'C, — (0,0) = Ci — (0, a), 


where (0, @) is a boundary identification vector proportional to the angular deficit. Here we 
subtract the identification vector (0,a) because this particular entwinement is counterclock- 
wise. For clockwise Wt we should rather add (0,+a). In that case points O, and Op are 
intersections of geodesics AC, and Cc. with faces of the static wedge. 

Applying the same procedure to the other entwinements, in a generic case we get a system 


of recurrent relations 
C = Ay(Ay'A ae (0, a)), (19) 
Co; = Aj(Aj*Co5_2 ot (0, a)), 
Con-1 = An(Ay B +- (0, a)), 


Co5-1 = Aj(Aj*Co544 + (0, a)). 
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Here A; = A; if the corresponding winding is W;. A; = Arr if the corresponding winding 





is W7,. In these formulae we pick up the upper sign if Wr, p and the lower sign if W; ;;. 
Note, that the Lorentz boost we have defined in terms of the embedding space coordinates 
acts non-linearly on the AdS3 points, therefore we can not simply expand the parentheses 
in (19). 
Then for each of the auxiliary arcs Co,C2441 we can derive coordinates of the refraction 
points Oo;, Oo,41, and write down corresponding lengths of the composing arcs. 


Later we will also show that not every formally generated sequence of windings does exist. 


IV. QUASIGEODESICS CONNECTING BOUNDARY POINTS 


To discuss causality properties of the dual boundary QFT, we will in particular need to 
consider boundary points with timelike separation. The conceptual problem we unavoid- 
ably encounter here is the absence of timelike geodesics connecting points on the conformal 
boundary of AdS. The equations for timelike geodesics can be derived from the following 
Lagrangian for a massive particle in AdS: 


2 





es 2) {2 sa 
(1+ R°)t Torte? il (20) 
Such a particle has two conserved momenta: 
E=(14+ Rt, (21) 
J = Ro. (22) 


Substituting them into the Lagrangian we obtain the radial equation of motion: 





&2 R? i 
“Te Tee ee ee) 
. 2 
R? =—(1+ R?) (+2) +e. (24) 


Clearly the right hand side of the equation turns negative as R — oo, and thus no real 
solution to this equation can exist. 

In the case of a stationary spacetime the obstacle could be easily surmounted by per- 
forming the analytic continuation of the metric to the Euclidean signature, calculating the 
Green function in terms of Euclidean lengths of the geodesics, and making the inverse Wick 


rotation back to real time. However in our case we deal with a spacetime that is not only 
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non-stationary, but wich has no good global notion of time. Hence we are forced to stick to 
the Lorentzian time. 

The way to implement the geodesic approximation for timelike separated boundary points 
in the single Poincaré patch has been suggested in [35]. Let’s turn for a second to the single 
patch of the AdS3 spacetime, covered by the Poincaré coordinates: 

dr? 
ds* = —r*dt? + —— ar dx". (25) 
r 
Again, a massive bulk particle has two kinetic invariants: 
Beart, (26) 


TS 1a, (27) 


but now we have two different classes of spacelike geodesics. 


e For J? > E?: 
r(A) = VJ? — E? coshr 
(A) = 29 + wp tanh (28) 
t(A) = to + x2 pr tanh A 

e For E? > J’: 


r(A) = VE? — J?sinhaA 
(A) = to — ge24r coth r (29) 
t(A) = to — pez cotha 


We will be interested in the geodesics of the second kind. As \ = 0 these geodesics approach 
the point r(0) = 0, which is the Poincaré horizon of the half-AdS chart. Regarding the 
horizon as a single infinitely far point (as in the theory of complex functions), we can consider 
two disconnected spacelike geodesics possessing the same kinetic invariants F and J, but 
emerging from two different timelike separated boundary points A(t4,x7,4) and B(tz, xz), 
as two branches of a single geodesic reaching the spatial infinity and returning back to the 
boundary. The length of such a geodesic will be divergent not only as r + oo (the standard 
holographic UV divergence), but also as r + 0, but this can be cured by an appropriate 


renormalization} The resulting expression for the renormalized length of the geodesic is 





3 For details see App. B and D of 
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simply 
£=In ((tg —ta)? — (ta —25)’) , (30) 
which gives the correct answer for the two-point correlation function of (1+ 1)-dimensional 


CFT 


1 
G(ta,ta;tp,tg) =e = ,§ #£P SA. on 
(ta A; lB B) ((tz — ta)? — (ta —2xp)2)> ( ) 


In the global AdS' space-time the Poincaré horizon has no special physical meaning, but 





we can still try to generalize this procedure to this case. 





The boundary field theory now is defined on S' x R! spacetime, and the two-point Green 
function that we must be able to reproduce via the geodesic approximation has the form 


: 


1 
Gta; Oa; te, ep) — | cos(tp = ta) = cos(bp = ba)|4 7 (34) 


Note that this function is periodic both in angle and time. While the angular periodicity 





is obvious by construction, periodicity in time emerges because of the finite size effects: an 
excitation created at some point in space and time starts dissipating, but later recollects 
and revives due to the spatial periodicity. 

We will need a function that defines angular separation between boundary points while 
properly maintaining the rotational invariance of the system. For instance, given two angular 


coordinates @2 = im and ¢ = 7, the difference between them along the shorter arc is 
D(a, 61) = —5 # ba ~ 1. (35) 
Thus we should use 
D(¢1,¢2) = mod (¢2 — d+ 7,27) —7. (36) 


Analogously, for the arithmetic average of two angular coordinates (that provides a point 


exactly at the middle of the shorter arc between ¢; and ¢9): 
1 1 
U(¢1, G2) = 5 (« + 2 — 276 (- COS (5 (oy — é2)))) (37) 


4 Tf Ax? > At?, the renormalized length is 





t= In (tp tay + (2A tp)’), (31) 


and the full Green function is 





G(ta,2a;te, 2B) = (32) 
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Note that the Green function has a symmetry: 


G(ta, 6a; tz, bp) = G(ta, ba; te +7, OB +7) ’ (38) 


where points on the r.h.s. can be spacelike separated while points on the |.h.s. have timelike 


separation: 
(tp —ta)? > D(dp,¢a)*?, but (tg —tat7)? < D(idp tm, da)’. (39) 


This symmetry can be used to construct a disjointed spacelike geodesic, with two branches 
reattached at the Poincaré horizon, connecting timelike separated points. 

If we represent the global Ad.S3 space-time as a cylinder, the Poincaré horizon consists of 
two planes cutting the cylinder at 45°. The orientation of the planes (as a rigid construction) 
can be chosen arbitrarily. Then consider a spacelike geodesic emerging from boundary point 
A (see Fig|6), and terminating at boundary point B*. Somewhere in the bulk it has a 
turning point P; where its radial coordinate R* = R(0) is minimal. 

Since we are free to choose the location of the Poincaré horizon, we can always orient it 
in such a way that the AB* geodesic intersects it at the turning point P,. This point P; can 
be identified with a point P, = P; + (m,7,0R), located on the other cutting plane. The arc 
P,B* can be then rotationally translated to this point: P, — P), and then B* + B, where 
B= B* +(a,7). Hereafter we will use disjointed “quasigeodesics” of this AP; P,B type to 
connect timelike separated points. 

The explicit analytic expression for the quasigeodesics can be derived in the following 
way. Let’s take the Poincaré chart geodesics (29), and rewrite them in coordinates of the 


embedding spacetime. The corresponding coordinate transformation is given by: 


r 


1 
Sg an); (40) 
A = ar, (41) 
2_ rl 2 42 
xX =o tre —t*), (42) 
=r, (43) 


Substituting here (29), and recalling the relations between the global and embedding coor- 
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FIG. 6. A spacelike quasigeodesic connecting two boundary timelike separated points A and B. A 


spacelike curve emerges from point A and reach the Poincaré horizon at point P,. Then it jumps 


to a mirror point P2 = P; + (m,74,0R) on the second plane of the horizon, and proceeds further 


to B. The arc PoB is a rotation of P, B* by az. 


dinates (3). we obtain 


—] EE? — 2 
V1+ R? cost = ( eesinh a, 


=> E 
1+ R?2 sin t = ———— cosh 
VE? - f 





ni 
EO eee 
_1_ Fp 2 
rt ee oie eee 





2V/ bE? — J? 


(44) 
(45) 
(46) 


(47) 


These can be solved to give us the embedding of the Poincaré chart spacelike geodesic into 
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the global AdS: 








2k 
t(A) = arctan [ee coth A) + to, (48) 
-1-E°+ / 
g(A) = arctan (+ tanh A) + do, (49) 
P 1S? . 3 
= ,/——~— h°A. 
R(A) |p > cosh* \ + ME? — 2) sinh” A (50) 


This quasigeodesic already has a 7-jump in time at the turning point \ = 0, and as explained 


before we also need to adjust the discontinuity in angle: 


~ (+ 


(A) = oA) + 70(A) = arctan ay tanh A) + 70(A) + do. (51) 


From now on we will omit the tilde. 


The kinetic invariants can then be expressed in terms of the boundary coordinates: 


+ to-ti 
sin 5) 








pe ; (52) 
sin Pea n.G) — cos eh 
D(¢2—n,¢1) 
cos Set 
es z (53) 
sin Diga—n1) — cos 254 
The integration constants tg and ¢ in can be represented as: 
1 
to = 5 (ti + ta) , bo = U(¢1,¢2 —T). (54) 


Inverting equation (50) we obtain dependence of the affine parameter on the radial coor- 


dinate: 





. A(E2 — J2)R2 — 4) 
= h 
A(R) = + arcsin ic EPP EAP? (55) 


where the minus sign is taken on the first branch of the geodesic or quasigeodesic (i.e. before 





the turning point, - when particle moves away from the boundary), and the plus sign is taken 
on the second branch (when particle moves towards the boundary). This function can be 


used to define the geodesic length, which is simply 


L£(R1, Ro) = Ax(R2) — Az (Fr), (56) 








for two points with radial coordinates R, and Ro. 


Note that for a geodesic connecting two boundary points the length is divergent: 








ee 4(B2 — J2)R2—4J2 
eee eA) 2 ic E+ P2+4PR eD) 


a 


and needs to be renormalized. The natural way to do it is to subtract the parameter 


independent divergent part, and define the geodesic length as 





(58) 


. 16(E? — J?) 
Ln im (A(R) — A(R) — 2In R) =In ( ) 


E4 — 2b7(-14 J?) + (14+ J?) 
The argument of the logarithm is always positive for quasigeodesics connecting timelike 
separated points, but can be less than 1. It means that the geodesic length after the 
renormalziation in principle can be negative. Using along with we obtain the 


correct result for the retarded Green function: 


1 
= hia. = 
Ge(ti, b1; t2, G2) = € ~ (cos(tg — t,) — cos(¢y — ¢1))4’ 


where (tg — t1)? > D(¢1, b2)", 





(59) 


where D(¢,, 2) is the function introduced in (36). The possible negativity of the renormal- 
ized geodesic length is the reason why singularities of the correlator can be captured in the 
geodesic approximation. 


Here we must pause for a second and stress Lorentz non-invariance of (59). We define 





SO(2,2) isometries of AdS3 in terms of the Lornetz boosts of the embedding R?”) space 
(6). If we take two boundary points A and B, and act on them with a bulk isometry 
transformation A of this type, we will observe that it does not preserve the Green’s function 
(59): 

G.(AA, AB) 4 G_(A, B). (60) 


It is the fundamental difference between holography of a Poincaré chart and holography of 
global AdS. In the first case the bulk isometries induce Lorentzian boosts on the boundary, 
so the Green function of a dual boundary field theory is a relativistic invariant object (33). 
In the second case the isometries rather act as conformal transformations leaving the Green 
function covariant, i.e. invariant up to some coordinate dependent scaling prefactors. 

In the holographic language this is encoded in the fact that the renormalized lengths 
connecting boundary points are dependent on the choice of the reference frame. Below when 
we consider the DeDeo-Gott time machine geometry, we should be especially careful about 
this, since the geodesics there are combinations of Lorentz invariant and non-invariant terms 
as, for example, in (17). The proper way to deal with it is explained in the first subsection 
of Sec. 
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In the next section we will analyze lensing of the quasigeodesics on the conical defects 
and calculate the Green function of the dual field theory in presence of the closed timelike 


curves in the bulk. 


V. THE TWO POINT GREEN’S FUNCTION 
A. The algorithm 


In Sec. [ITI] we have discussed the general idea of using the geodesic approximation to com- 
pute the boundary Green’s function for the DeDeo-Gott geometry. Now we will formulate 


an exact algorithm for that. 


e Introduce coordinate system on the unfolding of the double-cone space in such a way, 
that the physical (unremoved) part of the boundary consists of two stripes covered by 


coordinate intervals: 





t € (—o0, oo) , (61) 
n Ld an Ad an Ad nm Ad 
ge (-5-fae-f+Pee)u(s- Wk eae +t). 


e Fix two boundary points A and B. For simplicity we can choose A = (0, —}). 


e Fix the total number of windings N that a geodesic of interest undergoes on the way 
from A to B. In our simulations we will not go beyond N = 4, because the higher-order 


contributions to the Green’s function are highly suppressed. 


e For the given A, B, and the number N, generate all possible 4 - 3%! sets of the 
complementary points {C,,---C2y} corresponding to different sequences of windings 
{W,,--:Wwy}. The (quasi)geodesics then consist of N + 1 arcs AC), C2C3, ..., Can B, 


each of which is just a (quasi)geodesic curve in empty AdS3. 
e Impose that each of the “odd” complementary points C2;11 belongs to the causal future 
of the previous “even point”: Coj41 > Co. 


Let us elaborate on what the reason to do so is. The most clear question we can 
ask is whether causal propagation of a signal from the future to the past is possible. 


To define the dual retarded Green function in presence of the CTC in the bulk, we 
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should recall that evolution of a particle moving in the bulk of AdS can be split 
in two parts: “physical” continuous motion along a timelike or a spacelike geodesic, 
and “topological” time jumps caused by entwinement around the conical defects. In 
the holographic language geodesic branches C2;C2;1; correspond to the continuous 
evolution, and C34, — C42 identifications - to the time jumps. In absence of the 
closed timelike curves a signal could causally propagate from A to B if B belongs to 
the future light cone of A: B > A. A natural generalization of this prescription for the 
time machine case is to impose that this should hold true for all “physical” segments, 


1.€. Coi44 > Co Vi. 


For each of the causal quasigeodesics, solve for the intersection points {O1,---Oon}. 
The easiest way to do this is to transform for each winding back to the rest frame of 
the corresponding wedge. For example, if branch C9;C>;,; intersects first the trailing 
face of the 2nd wedge, and then the leading face of the 1st wedge, we perform a Lorentz 
transformation of the branch to the 2nd rest frame, then untwist the wedge by a simple 
coordinate transformation ¢’ = ¢ —t, such that angular location of the face remains 


still in these co-rotating coordinates, and solve the equation 
(PG o¢Coies (Aza) = coe (A2)) mod 27 = or, mod 27 . (62) 
Then we repeat the procedure in the 1-st rest frame: 


($05:Coi41 Avi+1) = b Oo; Coiet (A241)) mod QT = OL, mod QT : (63) 


Make sure that all these equations have real solutions (otherwise discard the geodesic). 


Make sure that if a branch is not expected to intersect other faces within the physical 
region of the space, it actually does not (fake intersections within the removed part of 
the unfolding are allowed). In other words, if an arc O2;O2; emerges from the face L, 


and terminates at the face 77;, it should not have intersections with L;; and T7. 


Calculate the lengths of all inner segments of the geodesic (O2;O2;+1). They are finite 


by construction and equal to 
LOuOstis = A241 — A2i : (64) 
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FIG. 7. A quasigeodesic with winding number N = 2 connecting boundary points A and B. 


e Renormalize the lengths of the boundary segments AO, and Oo) B as they are diver- 


gent: 
1 ren 
Lao, = 5g LAC, AG, (65) 
1 ren 
LogyB — 2 ConB AQN - (66) 


e Calculate the renormalized lengths of Lao, and Lo,,g in the original frame. As 
mentioned in the previous section, the renormalized lengths are not Lorentz-invariant. 
So, while we are free to constantly switch between different reference frames in order to 
calculate lengths of the finite inner segments O2;O2;,1, the renormalized lengths of the 
two boundary segments must be calculated in the original frame where we define the 
Green’s function. In our case it is the “centre-of-mass frame”, where the two conical 


defects are symmetrically boosted. 
e Finally calculate contribution of the geodesics to the Green function: 
GAB)=) 6°", (67) 
k 


where the index k runs over the set of geodesics that satisfy aforementioned conditions. 
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In accordance with the described algorithm we subsequently account for contributions 
from higher winding numbers starting with N = In other words, we formulate a kind of 
“perturbation theory” with the number of entwinements as a control parameter. 

One property of this series expansion must be comment on. Each geodesic contributes 
to the Green function exponentially: 


Ce (68) 


For higher windings the number of internal segments Co;Co;1, of the geodesic grows linearly 
in N, and so does its renormalized length L,.¢,. Therefore the corresponding contribution to 
the Green function is exponentially small’ On the other hand the total number of possible 


topologically different geodesics scales as 
nue’, (69) 


i.e. grows exponentially. Therefore in principle these two effects can compete and we can 
not say a priori that the higher order contributions to the Green function are suppressed, 
and the sum over entwinements is convergent. If not, this could mean that our setup is 
unstable and undergoes a Hagedorn like transition. 

However there are three different reasons for it not to happen. Firstly, by no means all of 
these 4-3%—! winding configurations satisfy the causality condition: Cy;,; + C2;. Secondly, 
even if the causality condition for the set of complementary points is satisfied, the geometric 
structure of the geodesics becomes more and more complicated as the number of windings 
increases, and it becomes hard to force a geodesic curve to undergo the concrete sequence of 
windings (it is easy to see on Fig|8{b)). Finally, the “decaying” exponent has a conformal 
dimension as a knob, so at large enough A it dominates over the “growing” exponent. 

Another way to understand convergence of the series expansion in all orders relies on a 
simple and general argument. Consider a germ of all possible quasigeodesics emerging from 
point A. The first segment of a generic quasigeodesic curve hits the boundary at some point 
C, first in the sequence of complementary points {C;}. If we go along the curve further, 
we will obviously see that it is defined uniquely up to the final moment when it reaches the 


physical part of the boundary at point B. Thus, for a fixed initial point A, for each of the 





> For negative times N = 1 windings do not contribute as they are due to the lensing on a single conical 


defect, that obviously can not lead to time travelling. But for positive times we take them into account. 


® Tn a generic case when Lren > 0. 
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FIG. 8. (a) An example of a physical N = 4 geodesic contributing to the Green function. (b) An 
example of physically impossible winding: for a given sequence of entwinements a geodesic can not 


be fit in the unremoved part of the spacetime. 


“first-in-the-sequence” complementary points C| the final point B is defined unambiguously. 
Now, as we have already emphasized, the singular contributions to the Green function come 
at the points where the renormalized geodesic length is infinitely negative, L’*” = —oo. It is 
possible if and only if C, is located exactly on a generatrix of the light cone emerging from 
point A, or C)N is located on the generatrix of the light cone of point B. Quasigeodesics 
having the complementary points C; and Coy right on the corresponding light cones form 
a zero measure subset among all possible quasigeodesics. Thus the set of boundary points 
where the Green function is infinite is also a zero measure subset of the boundary spacetime. 


Everywhere else the Green function is finite and well-defined. 


B. Phenomenology and discussion 


We are now ready to implement our computational algorithm for the DeDeo-Gott geom- 
etry. As shown in Sec. closed timelike curves in the spacetime are present when the total 


angular deficit is more than 27. For concreteness we impose a;7; = /3n, and the boost 
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FIG. 9. The Green function is sourced at point A with coordinates (0,—7/2). Blue lines depict the 
light cone emerging from this point. Any timelike line would cross the removed region and enter 
another strip (red dashed line). In order to avoid formulating the boundary field theory on both 
strips simultaneously we calculate the Green function on a timelike line very close to the light cone 


generatrix. 


rapidities Wy7; = +1. In their corresponding rest frames (in the co-rotating coordinates) 


the locations of the edges are taken to be 


a/2 bp, 
a/2+7 dnp, 


—a/2, 
—a/24+7. 


! 
Ih 


(70) 


/ 
L2 


It is more convenient to calculate the Green function also in the co-rotating coordinates: 
G"(t1, 15 ta, Oy) = G(ti, b1 — tj te, b2 — ta). (71) 


For simplicity we will mostly study the Green function on a one-dimensional timelike line 
passing through the point A. We should be careful here. Any timelike line originated in 
one physical strip crosses the cut out region and enters the second strip. To formulate a 
quantum field theory on both strips simultaneously is possible yet tricky due to the fact 
that on the unification of two parts of the boundary time can not be globally defined. To 


avoid this difficulty we will consider the Green function on a timelike line in a close vicinity 
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of the generatrix of the light cone. Then in a large range of times we will stay within one 
strip of the boundary. 


In other words, the object we will attempt to evaluate is (in the co-rotating frame) 
G" (0, —7/2;t, -7/2+ et), «<1. (72) 


We have performed the numerical calculation of the retarded Green function for negative 
times t < 0 up to N = 4 order, and for positive times t > 0 up to N = 2. 

Let’s discuss firstly the analytic behaviour of the Green function at negative times, - how 
the quantum particle behaves travelling back in time. A naive expectation would be to think 
that the Green function decays as t > —oo, and it is partially true. However as we can see 
at not very large negative times the function develops a number of non-trivial features, - 
peaks which we can interpret as the “most probable” regions of times the particle can reach 
using the time machine. The corresponding results are present on Fig{10} As an illustration 
we also provide a two dimensional plot for the leading N = 2 winding at negative times, 
Figll 

The origin of these peaks can be traced back to the fact that renormalized length of 
a geodesic can be negative. Generically at small conformal weights N = 2, N = 3 and 
N = 4 contributions are commensurate, but already at A = 2, higher entwinement terms 
are getting suppressed as compared to N = 2. However at specific points, where Leen < 0, 
the corresponding contributions to the Green function are getting enhanced in the large A 
limit!" forming a sharp peak. For instance, N = 2 set of geodesics contains such a curve 
around t2 = —1.9, and N = 4 set has a special point at t2 ~ —1.45, Fig/10} 

We have not performed numerical simulations for N > 4, but we can not exclude that 
such negative length curves can appear also at large N. The geodesic length is defined by 
lengths of internal segments (always positive) and lengths of the two boundary segments 


(that in principle can be negative): 


N 
TEM. -—. ren ren 
‘4B = L410, + LOswB + 5 LO 2005 * (73) 


i=1 


N 
If £45, + Lor pe < 0, and |LYG, + LOR Bl > do Loo:-202:-1, the Green function will get 


i=1 
a contribution that does not vanish in the large A limit. For a large number of internal 





” Strictly speaking, the geodesic approximation is reliable only in this limit. 
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FIG. 10. (a) N = 2, (b) N = 3 and (c) N = 4 contributions to the retarded Green function at 
negative times at A = 1.5. Discontinuities of the curves are artifacts of the geodesic approximation. 
(d) The retarded Green function at A = 1.5 (N = 2, N = 3 and N = 4 contributions are added up). 
For the large conformal dimensions peaks are enhanced, not suppresses, and we can see revival of 
the particle at moments preceding the excitation of the Green function. A not very large conformal 


dimension is chosen for convenience of presentation. Here ¢ = 0.1. 
segments it is not likely, but neither is impossible: while all internal lengths are finite, the 
renormalized negative lengths might be of an arbitrarily huge absolute value: 

40, <9, [L%6,| > 1, (74) 


thus dominating over positive contributions. 
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-3.5 -3.0 -25 


d - 2.0 -1.5 
FIG. 11. N = 2 contribution to the retarded Green function at A = 1.5 at negative times in two 
dimensions. The Green function is sourced at the red point (0,—7/2). We construct the theory 
only within one of the two strips of the boundary. Sudden break of the function signalizes that 


some regions of the spacetime are unattainable for the N = 2 quasigeodesics. 


In the case of a large conformal dimension it would mean that, if we were able to sum 
up contributions in all winding orders, the resulting Green function would have a shape of a 
comb with a number of peaks (in our calculations we discovered two of them). These peaks 
play a role of “pit stops” for a particle travelling in time, - they form a set of easily reachable 
coordinates in time. Hence we deal with specific “negative time” revivals. 

At positive times we have also discovered interesting features of the Green function. In 
the case of plain AdS3 geometry the dual light-like Green function (shifted away from the 
singularity) is decaying in time. In presence of the causality violating conical defects we 


detected a new peak of a high weight, Fig/12] signaling a revival of the excitation. 


VI. PHASES OF THE BOUNDARY FIELD THEORY 


We have calculated the Green function numerically up to N = 4 entwinements for the 
time machine geometry with a = 3m and w = 1. However it would be interesting to study 


how the properties of the Green function change upon changing the strength and rapidities 
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G(0, —7/2; t, (1 — e)t — 1/2) 
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1.0 1.5 2.0 


FIG. 12. The retarded Green function at conformal weight A = 1 at positive times. The plot 
demonstrates how the non-trivial N = 2 windings modify the original Green function (N = 1 
does not contribute when we consider the boundary theory within a single strip). The offset plot 
represents the N = 2 winding contribution separately. We have made the Green’s function timelike 


to avoid dealing with the light-cone singularity. Here « = 0.1. 


of the conical defects. We constructed the leading order N = 2 contribution to the retarded 
Green function at negative times for a € (1.17, ...1.957), and w € (0.1,...1.5) with stepping 
Aa = 0.057, Aw = 0.05., paying special attention to the location in time and and strength 
of the revival peak. 


The results can be schematically summarized in a form of a phase diagram, Fig]13} 


e If for a given value of a the rapidity w is not large enough to prevent the system from 


collapsing; the DeDeo-Gott geometry is forbidden (blue). 


e If for a given value of a the rapidity w allows for the existence of the DeDeo-Gott time 
machine, but still not very large, we clearly see the effect of revival, and the peak is 


sharper the closer 7 is to the lower bound (yellow). 


e If the rapidity is too large, the causality is violated, but excitations just decay and do 


not revive at negative times anymore (green). 
e At very small values of a the retarded Green function does not exhibit any non-trivial 
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FIG. 13. The four different “phases” of the boundary field theory. Blue spots depict the region 
of forbidden geometries, yellow spots are for the phase of negative time revivals, green spots 
form the region where the retarded Green function is non-zero at negative times, but does not 
exhibit reviving peaks in the large A limit. Red spots are where at the leading N = 2 order 
the boundary field theory retarded Green function does not demonstrate causality violation (i.e. 


G (0, —7/2;t, —1/2) =0, t < 0) despite the presence of the CTC in the bulk. Everything is based 





on the numerical simulations of the leading N = 2 contribution to the retarded Green function. 


We expect higher order corrections to change the diagram qualitatively, but not quantitavely. 


features at negative time even in presence of the closed timelike curves (red). However, 
this feature is likely just an artifact of N = 2 approximation, and we do not expect it 


to be there for higher windings. 


The profiles of the Green’s function at negative times are presented on Fig[14{a,b) for 


a = 1.57 and a = 1.77 respectively. The fact that revivals are seen only at not very large 
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FIG. 14. The negative time N = 2 contribution to the retarded Green’s function at (a) a = 1.5, 
w = 0.45, 0.50, 0.55, 0.60, and (b) a = 1.7, ~ = 0.75, 0.80, 0.85, 0.90. Here A = 1. On each of the 
plots the first three peaks are getting stronger and sharper as A —> oo (the yellow region on the 
phase diagram), while the fourth one is getting suppressed in the same limit (the green region), 
so we do not consider it to be an actual revival of a non-causal excitation. Counterintuitively, the 


weight of the Green function drastically decreases when we increase a or/and w. 


rapidities (and the effect is stronger as closer w to its minimal possible value) is surprising 
and contrasts to how causality is broken in the bulk. The stucture of CTC is defined by a 
and w, and the time jumps become stronger as the angle defects and rapidities are increased. 
Thus we rather should expect that for high ~ the time travelling along the CTC is more 
efficient in the sense that amplitudes of the classical free Green’s function defined on the 
boundary are getting enhanced as a w grows. In the interacting holographic dual field theory 
the retarded Green’s function is damped for larger a and w, so we can claim that causality 


in the boundary field theory is broken mildly as compared to the bulk. 


Another interesting feature of this system is that while the overall weight of the Green 
function drastically decreases when the rapidity w is taken away from the “forbidden region” 
on the diagram, the actual past time penetration depth (i.e. the deepest reachable point at 
negative times where G°"(0, —7/2;t, —7/2) #4 0) increases (though very moderately), and 


this is in agreement with the “naive” intuition. 
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VII. CONCLUSIONS 


In this paper we have analyzed properties of a two-point Green function in a (1+1)- 
dimensional field theory dual to the DeDeo-Gott time machine geometry. Using the geodesic 
approximation we have shown that AdS/C'F'T is capable of describing a quantum field theory 
when causality is violated, and have shown that the corresponding boundary propagator has 
remarkable features. We discovered that in presence of closed timelike curves in the AdS 
bulk a causal propagation of an excitation from the future to the past is possible on the 
boundary, and the retarded Green function exhibits peaks at certain negative times. At 
positive times analytic structure of the Green function also changes, and new singularities 
arise. 

Surprisingly, we have found that as we increase the strength of the conical defects a and 
the rapidity w, the causality violation in the dual field theory is getting milder in the sense 
that the weight of the retarded Green’s function at negative times decreases. 

Contra to the previous results on the dynamics of physical systems in time machine 
backgrounds [24 our calculations have demonstrated that sometimes it is possible to 
define evolution of an interacting theory in a time machine without imposing any additional 
self-consistency constraints. Despite the explicit non-causality the Green function does not 
have any uncontrollable pathologies. 

Our considerations leave a number of open questions. First of all, we have to understand 
how to interpret the boundary state dual to the DeDeo-Gott geometry, - whether this quan- 
tum state is pathological or just exotic yet physical state. From the boundary point of view 
a single conical defect, if its angular deficit is a = 27(1 — 1/N), can be thought of as state 
created by a non-local twist operator in a conformal field theory [84]. But what it means to 


have such an interplay of two independently boosted defects has to be clarified. 


Another thing we have not touched on in the paper is the entanglement structure of the 
boundary state. We focused on the properties of the retarded Green function, and thus 
analyzed the timelike quasigeodesics. However, even below the a = 7 threshold, when the 
CTC are not present in the system, due to the lensing it is possible to connect timelike 
separated boundary points just by standard continuous spacelike geodesics. If a certain 
generalization of the Ryu-Takayanagi conjecture is true in this case, it would mean 


the boundary state is timelike entangled [38]. Possible physical outcomes of this fact is an 
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interesting direction for future research. 
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